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We calculate the dc supercurrent through a Josephson tunnel junction consisting of an antifer- 
romagnetic (AF) or ferromagnetic (F) interlayer sandwiched between two d-wave superconductors 
(d). Such junctions exhibit a rich dependence of the Josephson current on the interlayer parame- 
ters, including the possibility of — tt transitions with varying temperature or interlayer thickness. 
Furthermore, we study d/l/d junctions when the dSC leads include subdominant magnetic corre- 
lations. Induced magnetism near the interface can strongly diminish the critical current for 110 
oriented junctions whereas no suppression exists for the 100 orientation. This may help resolve a 
long-standing puzzle of the critical current versus grain boundary angle in high-Tc superconductors. 

PACS numbers: 74.45.-(-c, 74.50.-fr, 75.50.Ee, 74.72.-h 



I. INTRODUCTION 

Interfaces and Josephson junctions between supercon- 
ductors and magnetic materials can generate novel low- 
energy spin-dependent Andreev bound states leading 
to highly unconventional quantum transport properties. 
For instance, junctions consisting of s-wave supercon- 
ductors (s) and ferromagnetic (F) metals have attracted 
great interest in recent yearsi. Such junctions have been 
shown to exhibit so-called — tt transitions^!^, where, de- 
pending on the temperature T and the width L of the 
interlayer, the ground state is characterized by an inter- 
nal phase difference of tt between the two superconduc- 
tors. This effective negative Josephson coupling is similar 
to what can happen when tunnelling through magnetic 
impurities^.. The possibility of — tt transitions may be 
utilized as a basis for future quantum qubits^, constitut- 
ing an important example in the field of superconducting 
spintronics^. 

Another promising situation involves interfaces be- 
tween antiferromagnets and superconductors. In this 
case spin dependent quasiparticle reflection at the AF 
surface, so-called Q-reflection, combined with Andreev 
reflection on the superconducting side, can lead to new 
low-energy bound states with important consequences 
for the associated proximity effect^i^. Furthermore, in 
s/AF/s Josephson junctions these bound states can en- 
hance the critical current Jc at low T^, and lead to 0- 
or TT-junction behavior depending on T and thickness of 
the AF interlayer^. For s/AF/s junctions the — tt be- 
havior is a true even-odd effect arising from qualitatively 
different spectra of the Andreev bound states caused by 
different symmetries of the odd and even AF interfaces^. 

We study the Josephson current through in-plane 
d/ AF/d tunnel junctions. Such junctions have not been 
studied before theoretically. Interestingly, our results 
are also relevant for d/F/d junctions. Based on both 
analytical calculations and numerical solutions of the 



Bogoliubov-de Gennes (BdG) equations, we determine 
the criteria for — 7r-junction behavior and predict un- 
usual T dependence of the critical current Jc{T). 

Intrinsic d/AF/d junctions may already be present in 
the case of high-Tc grain boundaries (GB) which induce 
AF surface states. Below, we also study the critical cur- 
rent through GB by modeling them as d/l/d junctions, 
where I is an insulating layer but where the leads contain 
subdominant magnetic correlations which become impor- 
tant near order parameter-suppressing interfaces. Both 
kinds of junctions mentioned above are cousins of the 
unconventional d/l/d junctions with uncorrelated leads 
which exhibit an unusual l/T behavior of Jc{T) at low 
T as well as possible (depends on misorientation angle) 
T- induced 0— tt transition o^°d^ . The experimental obser- 
vation of these effects is notoriously difficult due to the 
complexity of the barrier interface, characterized, in par- 
ticular, by facetting, twins and especially by many high 
transmission channels. Only recently have the main fea- 
tures associated with mid-gap state contribution to the 
Josephson current been observed in experiments^^'— i^. 



II. MODEL 

The Hamiltonian is defined on a two-dimensional (2D) 
square lattice (lattice constant a = 1) 

- ^nfiia +Ymi{hi^ - fiii) . (1) 

ia i 

Here, cj^ creates an electron of spin a on the site i, t 
is the hopping matrix element, jj, is the chemical poten- 
tial, and and rrii denote the superconducting and 
magnetic order parameters, respectively. The associated 
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BdG equations are given by 
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where ICfj — —tS^ij) + {ztarrii — n)5ij, with tr = +1/ — 1 
for up/down spin and 5ij and (5^^^^ are the Kronecker 
delta symbols connecting on-site and nearest neighbor 
sites, respectively. The net magnetization is Mi = 
rrii/U = i [(nil) — {nil)], ^^'^ ofF-diagonal block 2?^ 
describes d-wave pairing Vij = —Af^S^ij), where Afj = 
—V{ci^Cji — CiiCji) /2. The coupling constants U (V) are 
non-zero on (off) the L atomic chains constituting the AF 
interlayer. By Fourier transform parallel to the interface, 
we obtain an effective ID problem at the expense of in- 
troducing an additional parameter ky. The dc Josephson 
current jVr' between two neighboring sites r and r' is ob- 
tained from = ~{iet/h)Y,„ (cla^r'a) - {c^^Cra) ■ 

For more details on the numerical and analytical ap- 
proaches, we refer the reader to Refs. 



III. RESULTS 

For s/AF/s junctions, the — tt behavior as a func- 
tion of interlayer thickness L exists both for 100 and 
110 orientations^^. This is not the case for d/AF/d junc- 
tions, where the 100 case displays only 0-junction char- 
acteristics with an Ambegaokar-Baratoff like dependence 
of Jc{T). Therefore, we focus on the more interest- 
ing 110 oriented d/AF/d junctions. We discuss only 
identical (and identically oriented) junctions, and re- 
strict ourselves to the tunnelling limit where the current- 
phase relation is sinusoidal, and Jc = J{n/2). The 
110 oriented d/AF/d junctions are categorized further 
into d/AFeven/d and d/AFodd/d junctions depending on 
whether the interlayer consists of an even or odd num- 
ber of chains, respectively. In Fig. [IJa) and [IJb) we 
show typical self-consistent results for Jc as a function 
of T for even and odd interlayer-chains, respectively. As 
seen from Fig. [IK a), d/AFeven/d are 0-junctions with 
a 1/T-like dependence of Jc in the large-C/ limit. The 
small dip in Jc at low T is caused by the finite width of 
the interlayer and disappears in the limits ^/ L, U —i- oo. 
As shown in Fig. [TKb), Jc(T) in 110 d/AFodd/d junc- 
tions exhibits a surprisingly rich T dependence: as U is 
increased, the pure 0-junction at low U becomes a tt- 
j unction at high T, crossing over to 0-junction behavior 
at some T* which eventually vanishes in the large-C/ limit 
where Jc{T) ~ — 1/T. The systematic — 7r-junction os- 
cillations versus interlayer thickness L is shown in Fig. 
[IJc). The fcj,-resolved current corresponding to parame- 
ters similar to the green curve in Fig. [Ijb), is shown in 
Fig. [ijd). The same momentum region contributes to 
the current at all T , a fact which will ease the analyti- 
cal interpretation presented in the next section. Results 
qualitatively similar to those shown in Fig. [T]can be also 
obtained for thicker junctions with smaller values of U /t. 
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FIG. 1: (Color online) Jc{T) vs T for a range of U for (a) 
L = & and (b) L — 5. For clarity, the curves have been 
normalized to give the same absolute value at the maximum as 
the curve with the smallest U. (c) Jc at T = as a function of 
interlayer thickness L. The red squares correspond to negative 
Jc and have been multiplied by —1 in this semi-log plot, (d) 
the fcy-resolved current shown for L = 5 and U — lOt (along 
green curve in (b)). For all results shown here: jj. — and 
V — t. The latter leads to ^/L ~ 1 — 2 in agreement with the 
short in-plane coherence length of cuprate superconductors. 



We now study d/AF/d junctions within a quasiclassi- 
cal approach where, as usual, all characteristic energies 
are assumed to be much less than the Fermi energy Ep. 
We assume for the coherence length ^ :s> a,L, and the 
junction properties are conveniently expressed in terms 
of the scattering S matrix containing the reflection and 
transmission amplitudes. The Josephson current is car- 
ried entirely by phase-dependent Andreev bound states. 

For the 100 orientation the d-wave order parameter 
does not change its sign in specular reflection, but it 
changes sign in Q-reflection processes. By contrast, in 
the 110 case the d-wave order parameter changes its sign 
both in specular and in Q-reflection events. An impor- 
tant manifestation of this physical difference between ef- 
fects of Q-reflection for different AF-interface orienta- 
tions is that the — tt transition does not take place for 
100 orientation in d- wave junctions, but the transition is, 
in general, present in 110 d-wave, as well as 100 s-wave 
junctions. More formally, in the 110 case the specular 
and Q-reflection possess identical outgoing group veloci- 
ties and form the outgoing flow along one and the same 
direction. This permits the reduction of the problem to 
a standard situation with conventional number of incom- 
ing and outgoing waves, which determines the rank of 
the S-matrix. This is not the case for the 100 orienta- 
tion, when specular and Q reflection should be considered 
separately. This increases the rank of the S-matrix and 
makes ultimate results for 100 junctions with finite trans- 
parencies strongly different compared to the 110 case. In 
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the following we focus solely on the 110 orientated inter- 
faces. 

For d/AFodd/d junctions, the general structure of the 
S matrix is similar to that of d/F/d junctions with sym- 
metric F interfaces. This follows from the fact that in 
the (110) orientation and for an odd number of chains 
in the interlayer, all spins are aligned in the outermost 
chains. For (110) d/AFeven/d junctions, the outermost 
chains of the AF interface have opposite spin polariza- 
tions (but still all aligned within each chain) and the S 
matrix is isomorphic to the three-layer FIF interface with 
antiparallel orientations of the two F layera^^. The exis- 
tence of the even-odd (0 — tt) behavior shown in Fig. [ijc) 
follows directly from this link between the S matrices for 
110 d/AF/d and d/F/d junctions^. However, in order to 
understand the T dependence of Jc{T) and obtain quan- 
titative criteria for transitions between 0- or 7r-junction 
behavior, we turn now to the explicit calculations. 

Consider first (110) d/AFeven/d junctions where the 
transparency coefficients satisfy = Dg = D = 1 — R, 
resulting in the following Josephson current 



J(x,T) = ^^:^^tanh 



|A^|7 ' 
2T 
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where 7 = (i?sin^ ^ + D coa^ f ) ^^'^ X is the phase 
difference across the junction. Here, not only A'' and 
D, but also the spin- mixing parameter O (sinO(A:j,) — 

[m/2tcos(fcy/V2)] |l+ [m/4tcos(A;y/V2)]^| \ are 

all fcy-dependent, and the total current is a sum of Eq.Q 
over all values of ky-. However, as seen from Fig. [IJd), 
the fcj^-sum is unimportant for understanding the qual- 
itative behavior. Eq. ([3]) is valid for arbitrary trans- 
parency, and the critical current Jc{T) is plotted in Fig. 
HJleft) for a range of 8. In agreement with Fig. [IJa), 
the junction is always a 0-junction. Near Tc, Eq. ^ 
reduces to J{X;T) = e| A'^pD sinx/2r which coincides 
with the result in nonmagnetic (110) symmetric d/l/d 
junctions. However, at low T the current ([3|) is given by 
J{x,T) = e|A'^|Z)sinx/7 which, in the tunnelling limit, 
reduces to J{x,T) = e| A'^jZ? sinx/|sin ^ I- Therefore, 



due to the factor sin ■ 



we find the remarkable result 



that the current substantially exceeds the critical cur- 
rent in nonmagnetic (110) d/I/d junctions with the same 
transparency coefficient D (Fig. [2|. 

Next we discuss the T dependence of Jc{T) for 110 
d/AFodd/d junctions. As argued above, this junction 
is similar to 110 d/F/d junctions with symmetric F in- 
terlayer. In the tunnelling limit we obtain the following 
expression for Jc{T) (a = ±1) 
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FIG. 2: (Color online) Jc{T) obtained from quasiclassical cal- 
culations in the tunnelling limit D = 0.001 for d/AFeven/d 
(left) and d/AFodd/d (right), plotted for different spin- 
mixing angles (from = (blue) to O = tt (red) in intervals 
of O.OStt). The ky dependence of Q has been neglected here. 



which is plotted in Fig. [2l^right) for a = 1 (see be- 
low). The result for an arbitrary transparency can be 
obtained along similar lines to Ref. [l^. In the absence 
of magnetism, when = 0, a = — 1, and for zero trans- 
parency Z) — > 0, there are zero-energy Andreev bound 
states at both d/1 surfaces of the (110) d/l/d junction. 
With increasing 8, the mid-gap states on each d/1 sur- 
face evolve into spin-split Andreev states on a d/Yl sur- 
face. For a given ky, the energies of these spin-split states 
are Ed = ±A sin (I). This is different fr om the s-wave 
case where — ±A* cos (^)^, and therefore the behav- 
ior of the Josephson current in s/F/s tunnel junctions^ 
strongly differs from d-wave magnetic junctions. Eq.Q 
can be qualitatively understood as follows: in tunnel 
junctions the surface states ea further split and become 
phase-dependent due to a finite transparency. As a re- 
sult, four current-carrying interface Andreev states exist 
for a given ky. Eq. (jj]) represents the Josephson current 
carried by these states in the tunnelling limit, when two 
spin-split surface states on each side of the junction only 
slightly overlap through the interlayer. 

In the limit of a nonmagnetic interlayer (O = 0, a = 
— 1), only the second term in Eq.Q survives and one 
obtains Jc(T) = e| A''|2£)/2r, with the well-known 1/T- 
behavior for d/I/c? junctions. This result is the tunnelling 
limit of the more general current-phase relatio n^^i^^ 



J(X, T) = 2e I A''! sin | tanh 
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Hence, there are no — tt transitions in (110) d-wave non- 
magnetic junctions. This, however, is not the case, in the 
presence of magnetic interlayers with finite spin-mixing 
O. Finite values of 8 result in the appearance of the 
additional (first) term in Eq.((4]), which is comparatively 
small for small 8, and for d <-k /2 has the opposite sign 
compared to the second term. The second term in Eq.Q 
is, in its turn, strongly modified due to finite 8 at suf- 
ficiently low T . Indeed, it becomes exponentially small, 
if T is much less than the spin-split finite energies of the 
Andreev states Ed- At the same time, the second term 
in Eq.Q dominates the current at higher T, for exam- 
ple, near T^. For this reason the — tt transition arises 
in magnetic 110 d/AFodd/d tunnel junctions under the 
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condition 8 < 7r/2, as a result of the interplay of the two 
terms with opposite signs in Eq.Q. In principle, the 
change of sign of the total current in Q takes place with 
varying T for any small value of 8, but effects of finite 
transparency confine the conditions for the presence of a 
— TT transition to not too small values of 8. 

For deriving the conditions for the presence of the — 
TT transition in the tunnelling limit, it is convenient to 
consider two limiting cases of Eq.Q, one near Tc and 

another at low T. Under the condition 2^^™^\^) 
T < Tc, Eq.(|3]) reduces to the simple expression 



MT) = ^ cos 8 



(6) 



which is suppressed by the factor cos 8 compared to the 
corresponding nonmagnetic d/l/d junction. Eq. ^ is 
valid, in particular, near T^. Under the opposite condi- 
tion T < ^ sin (f ), Eq.® becomes 



Jc{T)^ae\M'WD„Ds 
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which is suppressed by the factor |sin (^) | compared to 
nonmagnetic d/l/d junction. Comparing signs of Eqs. 
(O and ([7]) , it is evident that the — tt transition takes 
place with varying T when cos 8 > 0, that is for 8 < -1. 
For a = \ (which is the case for d/AFodd/d junctions) 
and 8 < -I the 0-state is the ground state of the junction, 
whereas the 7r-state exists near in qualitative agree- 
ment with Fig. [Tfb). Note that — tt transitions in s/F/s 
junctions happens when the opposite inequality 8 > ^ 
is satisfiedi^. We stress that our results (|4]), ^ and ([7]) 
describe also the current in d/F/d junctions. 

Above we have extracted the form of the current 
characteristics in d/AF/d junctions via both numeri- 
cal BdG methods and an analytical quasiclassical ap- 
proach. An obvious question is how well these two meth- 
ods agree. To this end, in Fig. [3] we plot again the 
BdG results (normalized to 1.0) for the case U — IQt 
with L = 6 (left) and L — 5 (right), and show self- 
consistent quasiclassical fits to these curves. Here sin 8 = 
m/2tcosky/\/2/[l + (m/4t cosfcy/\/2)^] and the special 
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FIG. 3: (Color online) Critical current Jc{T) vs temperature 
for (left) d/AFeven/d and (right) d/AFodd/d junctions. The 
green dots display the same U = lOt BdG data points shown 
in Fig. [T^,b (normalized to 1.0). The solid curves are qua- 
siclassical fits where the special choice cos ky/V^ = 0.7 has 
been taken, and the transparency D has been adjusted to fit 
the BdG data. 



choice cosky/^/2 = 0.7 has been taken. The trans- 
parency D has been adjusted to fit the BdG data. As 
seen, there is overall very good agreement. At low T some 
discrepancy can be detected, which we believe originates 
from the finite interlayer thickness used in the BdG cal- 
culations and/or the different bandstructures (circular vs 
square Fermi surface in the quasiclassical and BdG ap- 
proach, respectively). Disregarding any explicit ky de- 
pendence of the transparency coefRcients and the 8 pa- 
rameter in the quasiclassical calculations may also play 
a role. 

Experimental detection of — tt transitions in d/AF/d 
junctions may be possible in sandwich structures of high- 
doped and un-doped high- Tc materials similar to what 
was constructed for c-axis junctions^". Recently, Oh et 
al?^ , invented a spatially controlled doping method and 
fabricated in-plane 100 d/AF/d junctions. Our results 
show that the fabrication of similar 110 junctions and a 
detailed study of their Josephson current characteristics 
holds the promise of several unusual properties as shown 
in Figs. [l][21 Realistic junctions will contain regions with 
varying interlayer thickness, but if these are sufficiently 
few, the regions with shortest thickness should dominate 
the current. Alternatively, one needs to average the cur- 
rent over interface imperfections. Jc in even junctions 
dominates at low T only in the limit of large U. Then 
we have a small 8 and 0-junction with a low-T anomaly 
in Jc- Otherwise critical currents in even and odd junc- 
tions are of the same order. For 8 > 7r/2 (i.e. m < At) 
the currents have identical signs at all T (0-j unctions). 
For 8 < 7r/2, the 7r-junction state arises in odd junctions 
near resulting in an overall cancellation of odd and 
even contributions to the current. 



IV. GRAIN BOUNDARY JUNCTIONS 

Finally, we turn to the question of Jc through grain 
boundaries, where a strong discrepancy between theory 
and experiment has been known for some time: when 
the GB is modeled as a d/l/d junction the zero-energy 
state existing in the 110 orientation results in a large low 
T increase of Jc as compared to the 100 orientation (see 
dashed lines in Fig. d]). However, the opposite behavior is 
obtained in experiments: Jc is largest for 100 orientations 
and drops exponentially with increased angle between 
the GB and the crystal axis22,. We model the GB using 
Eq.([T]) in a d/l/d geometry with a potential V{ni^ + riii) 
inside the insulating layer (I) and [/ 7^ in the leads 
only. For sufficiently small U, magnetization is absent in 
the superconducting leads, but the magnetic correlations 
can lead to instabilities near interfaces that suppress the 
dSC order paramete r^'^i^^ , as shown in the inset of Fig. 
m The main body of Fig. d] shows Jc{T) for a range 
of U all small enough not to induce magnetization in the 
bulk. Contrary to the 100 orientation, Jc through 110 GB 
can be significantly reduced by surface-induced magnetic 
order for T < Tm, where Tm is the critical temperature 



5 



„x10" 

0| — 1 — 



0.2 




FIG. 4: (Color online). Jc{T) versus T for 110 and 100 d/l/d 
junctions with AF correlations in the d leads which have dop- 
ing level X — 0.1. The I region is modeled with L = 2 and 
a potential of V = 30t. In the 100 case, the same curve 
(black circles) is obtained for all U < 2.1t. (inset) Example 
of surface-induced magnetization (blue) and the suppression 
of A** (red) at a 110 interface, shown here for U = 2.0t. 



ness and temperature using both numerical BdG diag- 
onalization and analytical quasiclassical methods. For 
an odd (even) number of antiferromagnetic chains in 
the interlayer, the current characteristics of 110 ori- 
ented interfaces display 7r(0)-junction behavior. In ad- 
dition d/AFodd/d junctions can exhibit vr — transi- 
tions as a function of temperature. We have shown that 
in terms of the spin-mixing parameter 0, the condition 
for the latter is given by < This is the opposite 
regime as compared to leads vifith s-wave pairing sym- 
metry where temperature-induced tt — transitions take 
place for 8 > ^. Another important difference between 
s/AF/s and d/AF/d junctions exists for the 100 orien- 
tation, where d-wave junctions always are 0-junctions 
whereas this is not the case for s-wave superconductors. 
Finally we studied grain boundary junctions modeled as 
d/I/d junctions but with subdominant magnetic correla- 
tions in the superconducting leads allowing for interface- 
induced magnetism near grains which tend to suppress 
the d-wave order parameter. We showed that this mech- 
anism can lead to larger critical currents for the 100 ori- 
entation than for 110, in qualitative agreement with ex- 
periments. 



for the surface magnetization. In fact, as seen from Fig. 
|4] there exists a range of U where Jc at low T becomes 
smaller in the 110 orientation compared to the 100. This 
shows the importance of competing surface effects even 
though a complete understanding of the physics of GB 
junctions requires more detailed microscopic calculations. 

V. CONCLUSIONS 

We have studied the dc Josephson current through 
d/AF/d tunnel junctions as a function of interlayer thick- 
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